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Abstract 

We investigate the effect of the stochastic gravitational wave (GW) background produced by 
kinks on infinite cosmic strings, whose spectrum was derived in our previous work, on the B-modc 
power spectrum of the cosmic microwave background (CMB) anisotropy. We find that the B-mode 
polarization due to kinks is comparable to that induced by the motion of the string network and 
hence the contribution of GWs from kinks is important for estimating the B-mode power spectrum 
originating from cosmic strings. If the tension of cosmic strings // is large enough i.e., Gfi > 10~ 8 , 
B-mode polarization induced by cosmic strings can be detected by future CMB experiments. 
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I. INTRODUCTION 



Cosmic (super) strings can be produced in the early Universe at the phase transition 
associated with spontaneous symmetry breaking Jl], the end of supersymmetric hybrid in- 
flation 2|, yj , or the end of the brane inflation |4|, |5| . They can be a clue to particle physics 
beyond the standard model and the history of the early Universe, which is difficult to obtain 
in terrestrial experiments. How to find signatures of cosmic strings in cosmic microwave 
background (CMB) experiments has been extensively discussed for decades. Especially B- 
mode polarization of the CMB induced by the cosmic string network was investigated in 



many papers HlO|. 

B-mode polarization, which has not been detected yet, is polarization of the parity-odd 
type. It cannot be produced by the primordial scalar perturbation from the inflationary 
era, which is widely believed to be the main origin of the present structure of the Universe. 
On the other hand, the tensor perturbation can be a source of B-mode polarization. Some 
inflation models can produce the intense tensor perturbation enough to generate detectable 
B-mode, while others cannot. 

Cosmic strings can also induce B-mode. Cosmic strings move in the Universe in a very 
complicated and nonlinear way, constantly generating all types of perturbations, scalar, 
vector and tensor ones. Therefore, dynamics of the cosmic string network induces B-mode 
and it reaches an observable level if the tension of cosmic strings, fi, is large enough, say, 
G[i > 10 -7 10J. Here, G denotes the Newton constant. 

In this paper, we point out that there is an additional source of B-mode when the cosmic 
string network exists. It is the stochastic gravitational wave (GW) from kinks on infinite 
strings. 1 In the previous paper ll|, we investigated GWs emitted from kinks on infinitely 
long strings, and found that GWs with a wavelength comparable to the Hubble horizon 
scale are generated. These long wavelength GWs can produce an observable B-mode in the 
CMB. As a result, we find that the contribution of such a GW background to a B-mode is 



1 The effects of the kinks on infinite strings are partially reflected in the calculation in Ref. [9[ based on 
the lattice simulation. However, such a simulation covers only the limited period of the evolution of the 
string network, so the correct kink distribution on infinite strings cannot be taken into account by this 
method. Moreover, it is impossible to completely separate the GWs from kinks from those emitted at 
the phase transition. Therefore, our calculation based on the kink distribution derived analytically is 
complementary to the calculation in Ref. 
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comparable to that due to dynamics of the cosmic string network and detectable by future 
CMB experiments, such as PLANCK or CMBpol. 

This paper is organized as follows. In section 2, we present the formalism which we 
adopt in order to compute the BB power spectrum. In section 3, we briefly review the 
result of Ref. ll| for the GW spectrum from kinks which will be used in the following 



analysis. In section 4, we show the resulting BB power spectrum and discuss its observational 
implications. Section 5 is devoted to summary. 



II. FORMALISM FOR COMPUTATION OF THE BB POWER SPECTRUM 



In this paper, we adopt the formalism described in Ref. 12|. The tensor mode of the 
metric perturbation is defined as gij(t, x) = a(t) 2 (5ij + £)y(i, x)), where a(t) is the scale 



factor, and it is symmetric, transverse and traceless: 
expand Dij in the following form, 



0. We 



x) = / dV q - x Ai(*,q)= I rfV q %(g,A)P(t,q,A), 

J A=±2 J 



(1) 



where q = q/|q|, A denotes the helicity of the GW and e^(g, A) is the polarization tensor. We 
can think of T>(t, q, A) as a stochastic variable since it is the product of random GW emission 
by kinks on infinite strings. Its root mean square is inferred from previous paper Uj and 
given in the next section. 

As described in 12j, this metric perturbation relates to the polarization of photons. We 



do not explain the detail here, but write down only several important equations. The BB 



13j, where aB,e m is 



power spectrum, C B b,£, is defined by (a B/rn a B ,e> m >} = Cbb^m'^. 
some combination of the coefficients of the multipole expansion of the Stokes parameters, 
Q{n) and U(h), by the spin-weighted harmonics, as/m can be written as 

,t 



Q>B,tm 



A=±2 J ,Jo 



A) 



>>P + P Z 



d \ je(p) 



dp J p 2 



(2) 



p=qr(t) 



where T is the present CMB temperature, to is the age of the Universe, r(t) = J t ° 

is the spin-£ unitary representation of the rotation group, S(q) is the rotation which 
takes the three-axis into the direction q and je is the £-th spherical Bessel function. P(t) = 



u c (t) exp(— J uj c (t')dt') is the so-called visibility function, which has sharp peaks at the 



moment of recombination and reionization. u c (t) is the rate of Thomson scattering. \1/ is 
the function, which satisfies 



vD(t, q ,A) = - fdt'e-ti dt ""^ 
2 Jo 



-2T>(t',q,\)K[q 



dt" 



t , a(t") 



u c (t')F q 



dt" 



K{x) = j 2 (x)/x 2 , F(x) = j (x) - 2j 1 (x)/x + 2j 2 {x)/x 2 



*(*',q,A) 
(3) 
(4) 



The definitions of the Stokes parameters, as/m and \1/ are found in 12j. If we know the way 
for T> to evolve precisely, we can get \1/ through Eq. and calculate the power spectrum 
by integrating Eq. ([2]). However, we cannot know the phase of V, which varies randomly, 
since the stochastic background of GW is formed by random and continuous accumulation of 
GWs from kinks. We can find only the expectation value of its amplitude. Nevertheless, we 
can estimate the BB power spectrum using the 5-function-like property of P{t), as described 
in the Appendix. The B-mode power spectrum is calculated as 



C bb/ ~tt 2 T 2 / dqq 2 Mq)V 2 (t rcc ,q) + B e (q)V 2 (t mil q) 



(5) 



where T> is defined in the next section, A^{q) and B^{q) are defined in Appendix, and t re c(^rei) 
is the cosmic time at recombination(reionization), where the visibility function has a peak. 



III. SPECTRUM OF THE STOCHASTIC GRAVITATIONAL WAVE BACK- 
GROUND 



The amplitude of the tensor perturbation is found as below. In ll|, we derived the 
spectrum of the stochastic GW background, using this kink distribution function [j^j], which 
describes the abundance of kinks for a given sharpness. In the matter-dominated(MD) era, 
the energy density of GWs of frequency ~ u is 



dp 
dlnu 



lOGfi 2 (ut) c ™t- 2 for t- x < cu < J^ D \t) = (!f) Am r 1 

10G/i 2 rf) {uotfn- 2 for c! MD) (t) < u < J 2 MB \t) = (±f) Am (£) V 1 

(6) 

where A m = -0.8, A r = -0.92, C m = -0.17, C r = 0.14, D = 0.11, t eq is the cosmic time of 
the matter-radiation equality and t r is the time when the reheating completes. GWs with 
frequency larger than a4 MD ^(£) are irrelevant because their wavelength is too short to affect 



the large-scale density perturbation probed by CMB observations. The periods concerning 
the CMB polarization are only those around the recombination and the reionization, hence 
it is sufficient to consider the matter-dominated era only. This GW background consists 
of GWs emitted toward random directions from random points in the Universe at random 
time. Therefore, we can think of it as being isotropic and homogeneous. 2 

We want to connect this expression of the energy density to the amplitude of GWs T>. 
At the small scale where the cosmic expansion can be neglected, the energy density of GWs 
can be written as 



16] 



P = 7^(A J (i,x)Z^(t,x)>. (7) 



1 

32kG' 

This expression applies to GWs from infinite strings, whose wavelength is shorter than the 
Hubble radius. Under the present notation, the energy density of GWs whose frequency 
~ uj can be written as 

1 ;u 3 a 3 ti 2 (t,q). (8) 



dlnu; 2G 

Here, we set (D(t, q, \)V(t, q', A')) = V 2 (t, q)5 X \'5 3 (q - q') and used the fact that V(t, q) is 
oscillating with frequency q/a — u. Eventually, from Eqs. ([6]) and ([8]) we obtain 

VWGfiq-^ioutf™/ 2 ^ 1 for r 1 < u < J^ D \t) 

V{t,q) ~ { \ —D/A r . (9) 

^OGfiq^ 2 r-f) (utfr/H- 1 for W i MD) (t)<u;<u4 MD) (t) 

IV. BB POWER SPECTRUM 



Now let us calculate the BB power spectrum. As 



or the cosmological parameters, we 
Besides, we have to specify the 



used the result of the 7-year WMAP observation fl7 |. 
ionization history, or the shape of u c (t). It is given as u c (t) = <jTn e (t), where or is the 
Thomson scattering cross section and n e is the number density of electrons. We calculate 
the time evolution of n e (t) around the recombination epoch by using the RECFAST code jisl . 
Concerning the reionization, we make an approximation that the reionization occurs 



2 In Rcf. [11], we omitted GWs which do not overlap others from the "background", following the pre- 



scription given Rcf. [15|. Here, however, we do not consider this subtlety and include all GWs in the 
background, because we are paying attention to only long wavelength modes, most of which overlap 
others. 
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suddenly at some redshift z re : n e jumps from to some value n e o at z re . Thereafter, n e 
decreases in proportion to (1 + z)~ 3 . In short, we assume 

for z > z re 

u c (t) = { , (10) 

u c o(l + z)~ 3 for z < z ie 

rich is 



3. 



around the reionization. We set z Te = 10.4 according to Ref. [17] . u c o is a constant w 
determined so that f.° dtujJt) conforms to the reionization optical depth r = 0.087 

We show the resulting BB power spectrum in Fig. [T] (thick red) with the spectra produced 
by the string network dynamics (blue), the inflationary tensor perturbation with tensor-to- 
scalar ratio of 0.1 and 0.01 (black) and the lensing effect (green). We also show the sensitivity 
curves of PLANCK and two different realizations of planned CMBpol satellites, EPIC-LC 
and EPIC-2m. The spectrum induced by the string network dynamics is drawn by CMBACT 
20, [21] and those originating from the inflation are obtained by the CAMB code B. The 
value of G[i is set to be 10~ 7 , which is close to the present observational upper bound 23H25]. 



In computation using CMBACT, the network parameters are set as follows : the wiggliness 
a r = 1.8, the r.m.s. string velocity v r = 0.64, and the ratio of the correlation length to the 
cosmic time 7 r = 0.3, where subscript r means the values in the radiation era. They are 
derived from the results of the simulations 
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291 ]. and extrapolated to the matter era by 



the procedure in the code of CMBACT. 

We can see the spectrum induced by the GWs from kinks on infinite strings has two peaks, 
one of which is located at t ~ 100 and the other at I ~ 5. The peak at i ~ 100(£ ~ 5) 
is induced by GWs which exist at the recombination(reionization). At every moment, the 
lower limit of frequency of existing GWs is roughly the Hubble parameter at that time. 
Besides, the amplitude of GWs declines toward higher frequency. As a result, the position 
of each peak is set by the Hubble parameter at the recombination or the reionization. 
Remembering that GWs of frequency comparable to the Hubble parameter are, as discussed 
in 11J, emitted by new kinks, one finds that the peak at i ~ 100(£ ~ 5) is due to GWs 
emitted slightly before the recombination(reionization) by kinks which are born a little 
before the recombination(reionization). Since the BB spectrum by kinks is comparable to 
or somewhat grater than that by network dynamics in some regions, the total BB power 
spectrum induced by cosmic strings is deformed by the effect of the GWs from kinks on 
infinite strings. It is natural that the effect on B-mode of GWs from kinks and that of GWs 
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by global dynamics of strings are comparable. The network dynamics can induce B-mode 
also through the vector perturbation it produces . We expect that its contribution is also 
comparable to that of GWs from kinks, although strict comparison is difficult and requires 
numerical calculation. The magnitude of the spectrum is proportional to (Gfi) 2 , and hence 
when we take a different string tension the shape is unchanged but the whole spectrum 
moves upward or downward. If Gfi > (a few) x 10~ 7 , this spectrum can be observed by 
PLANCK, and if Gfi > 1(T 8 , it can be detected by CMBpol. 

The peak around £ ~ 5 associated with the reionization has a characteristic shape. 
However, it might be an artifact of the approximation that we put ^ out of the time 
integral in (T5]), assuming that P(t) has a sharp peak around the reionization 3 (see ( 1A1I) ). 
The actual peak might be smoother. In fact, B-mode is induced over the finite time around 
the reionization, and those which are produced at different moments have their peak at 
different i. It is expected that the total BB spectrum, which is the envelope of such peaks, 
has a smoother shape. On the other hand, we expect that FIG. 1 shows the correct position 

3 This approximation is much more valid around the recombination than the reionization, therefore the 
peak associated with the recombination does not have bump-like feature. 
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and height of the peak around l ~ 5. 



V. CONCLUSION 

In this paper, we have studied the effect of the stochastic GW background induced by 
kinks on the infinite cosmic strings on the BB power spectrum of CMB polarization. Using 



the GW background obtained in Ref. ll|, we have estimated the resulting BB power spec- 
trum. We found that this effect is comparable to that of the vector and tensor modes induced 
by motion of the cosmic string network and may leave observable signatures in the spectrum. 
If the cosmic string tension is large enough, the BB power spectrum by cosmic strings will 
be detected by future/on-going satellite experiments such as PLANCK and CMBpol. If it is 
discovered by the CMB experiments, then the direct detection of GWs from cosmic strings 
by pulsar timing arrays or space-laser interferometers may further confirm the existence of 



the cosmic string 



Appendix 

In this Appendix, we derive Eq. ([5]) using the 5-function-like property of the visibility 
function P(t). \I/(i,q, A) varies more slowly than P(t) at its peaks. Therefore, the time 
integral in Eq. ((2]) is approximated as 

° dtP(t)^f(t, q, X)xe(qr(t)) ~ ^(t Iec , q, A) / dtP(t) X e(qr(t))+*(t ieh q : \) [ dtP(t) X e(qr(t)), 



(Al) 

where / rec (/ rei ) represents integration around the recombination(reionization), and Xe( x ) = 
(8x + x 2 d/dx)(ji(x)/x 2 ). Then let us estimate ^{t rcc , q, A) and \I/(t re i, q, A) from 
Eq. ([3]). First, we consider \l/(t rcc , q, A). The factor exp(— f* rcc dt"uj c {t")) has the 
property that it rapidly increases from when t' approaches t rec . The combination 
u) c (t') exp(— J^ ICC dt"u c {t")) also has such a property, while T> and \1/ vary more slowly than 
these functions. Then, we obtain 

dt" \ 
, a{t") J 

^( trcc , q , A) jf dt'u c {t') exp (- dt"u c {t")^ F (a jf" . 

(A2) 



(l-trec \ / j'ticc 

- dt"u c (t")j K (q j 
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In contrast, exp(— J^, im dt"u c {t")) is almost 1 around t' = t re i, since u c does not increase 
enough around the reionization epoch. However, we can put T> out of the time integral 
because of the property of K J* t * rei ^pry j • K(x) decreases proportional to x~ 3 for large 
x. Besides, for q > a(t Tci )/t Tci , for which V(t rei , q, A) has a nonzero value, q J^ ro1 



{t") 

^aiy ^1 - (j-^j / J grows rapidly when t' goes away from t rei . After all, K (^q J^ C1 ^pr) 
has a sharp peak at t' = t re ;. Thus we can estimate ^(t Te i, q, A) as above, 

(t rci , q, A) ~ -3V(t ieU q, A) jf" dt'exp jf " d*"w c (*")) K L jf " 

+ ^(t rci ,q,A) / ^(tOexpf- [™ dt"u c (t"))F(q ' 



rei \ Jt' / \ Jt 



a(t") 



(A3) 

Connecting the above estimations, we finally get 

Cbb/ ~ vr 2 T 2 y dgg 2 [A £ (q)V 2 (t vcc , q) + B e (q)V 2 (t Ich g)] , (A4) 



where 



*J / /.#\ _ P*rec Ji.ll,. f+f'\ _ / 



2 7rec \ A' ^(t") 

and -B^(g) is the function which we can get by substituting t rec in Ag(q) for £ rei . 
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